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PERMANENT MAGNET AND SHIM DESIGN 

FIELD OF THE INVENTION 
This invention relates to the design of permanent magnetic structures for 
5 generating a specific magnetic field in a Region Of Interest (ROI) or modifying an 
existing magnetic field in an ROI using shims. 

BACKGROUND OF THE INVENTION 
Permanent magnets are of great value in many fields of science and 
engineering. In particular, they find application in the field of Magnetic Resonance 
10 Imaging (MRI) and Spectroscopy (MRS). Within these latter applications it is of 

paramount importance that extremely high degrees of magnetic field homogeneity 
(of the order of one part per million) are achieved within an ROL 
In permanent magnets the distribution of the magnetic material determines 
how the resulting magnetic field varies over space. In both permanent magnetic and 
15 electromagnetic MRI and MRS machines, manufacturing imperfections result in 
unacceptable inhomogeneities. To improve the homogeneity of an existing field 
over an ROI it is known to place permanent magnetic materials (termed passive 
shims) in the magnetic field to achieve this improvement of homogeneity. 

The problem of distributing magnetic material to produce, or modify an 
20 existing magnetic field is one of a class known collectively as the near field 
problem. In particular, it requires the determination of the distribution of material 
to produce a field of particular character - varying in a specific manner in space, for 
example a uniform field in a proximate region. The total magnetic field from a 
magnetic structure is a nonlinear convolution of the spatial distribution of material 
25 with the field of a magnetic dipole. While it is relatively straightforward to 
determine the resulting field from a given distribution of magnetic material, 
generally, the components of a convolution are not invertable. Consequently, 
determining the optimum distribution of magnetic material to produce a magnetic 
field of a specific type in an ROI has proven to be extremely difficult. 
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One known method of electromagnetic coil design is the Target field approach 
as disclosed in European Patent No. EP-A- 0 252 634 or its US equivalent No. 
4896129. As it is based on an assumed a priori mathematical expression for the 
total field, an optimum resulting field is hard to achieve. 

An alternative approach is disclosed in European Patent EP-A-0601 101 and its 
US equivalent No. 5266913, the content of which is hereby incorporated by 
reference. It is concerned with the design of shielded electromagnetic coils, and 
employs forward substitution and least squares minimization to optimize the field 
produced by electromagnets specifically within the ROI, whilst achieving a null 
field outside of the confines of the coil structure. 

There are significant differences in the form of the magnetic fields generated 
by active electromagnets, and those resulting from the use of magnetic materials. 
Consequently, electromagnetic coil design and the design of permanent magnets 
involve quite separate design considerations. 

The object of the present invention is to achieve an improved method of 
designing permanent magnetic structures. 

SUMMARY OF THE INVENTION 
In accordance with the present invention there is provided a method of 
designing a permanent magnetic structure for generating a desired magnetic field in a 
region of interest (ROI) comprising: 

-setting the dimensions of the permanent magnetic structure 
-setting a hypothetical distribution of magnetic material over the surface of 
said structure for each member of a finite set of spatial frequency harmonics of unit 
magnetisation amplitude 

♦calculating the magnetic field in the ROI from each member of said set 
-calculating amplitudes of said harmonics by the method of least squares in 
order to achieve the optimum approximation to the desired field in the ROI 

-scaling the members of said set according to said optimised harmonic 
amplitudes 
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-establishing, by summing said scaled members, the required distribution of 
magnetic material in order to generate the desired magnetic field. 

The permanent magnetic structure may either be a permanent magnet 
assembly for producing a magnetic field ab initio, or may be a permanent magnetic 
5 structure for modifying an existing magnetic field, by means of shimming. 

The method may also be applied iteratively in order to take into account the 
finite thickness of the required distribution magnetic material, which will inevitably 
have influence on the field. 

Further preferred embodiments of the present invention are presented in the 
1 0 dependent claims. 

BRIEF DESCRIPTION OF THE DIAGRAMS 
Figure 1 shows a schematic representation two parallel, finite, flat plates upon 
which magnetic material is distributed. 

Figure 2 shows an alternative distribution composed of a hollow elliptical, 
15 cylindrical structure upon which magnetic material is distributed. 

Figure 3 shows the magnetic dipolar field (Bd) produced by magnetic material, 
resolved into orthogonal components Bd r , Bd z and Bd G . 

Figure 4 shows schematically a hollow, cylindrical structure with a circular cross- 
section, upon which magnetic material is distributed. 
20 Figure 5 shows the permissible, even, cosinusoidal axial harmonics of the magnetic 
material distribution on a structure of finite length. 

Figure 6 is a plot of the axial magnetic field Br 2 against distance z resulting from an 
infinitesimal circular ring of uniformly magnetic material of unit magnetization 
calculated directly by numerical integration. 
25 Figure 7 shows a plot of the Fourier Transform (FT) in the axial direction (z) of the 
axial magnetic field of Figure 6. 

Figure 8 shows a plot of the FT, in the axial direction, of the axial magnetic field 
resulting from an infinitesimal circular ring of uniformly magnetic material of unit 
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magnetization evaluated by an alternative, analytical method, directly in the Fourier 
Kernel domain (k-space). 

Figure 9 shows a plot of the axial distribution of magnetic material that will achieve 
an optimal, uniform field, designed in accordance with a first embodiment of the 
5 present invention. 

Figure 10 shows a schematic representation of a magnetic structure designed in 
accordance with the invention, to generate the distribution shown in Figure 9. 
Figure 1 1 shows a plot of the resulting axial magnetic field in the axial direction 
produced by the distribution of magnetic material shown in Figure 9. 
10 Figure 12 shows a plot of the deviation of the above magnetic field from the 
desired, uniform one, in parts per million (ppm), over the ROI. 

Figure 13 shows a plot of the deviation of the above field over, and proximate to, 
the ROI. 

Figure 14 shows how a distribution of paramagnetic and diamagnetic materials on a 
15 surface can be arranged to produce a sinusoidal variation of magnetization along 
that surface. 

Figure 15 shows a plot of the amplitude of the optimum axial distribution of 
magnetic material over a cylindrical structure to provide a magnetic field that varies 
linearly in a transverse direction in the ROI, designed in accordance with the present 
20 invention. 

Figure 16 shows a schematic representation of how a shim magnetic structure 
designed in accordance with the invention and using diamagnetic and paramagnetic 
materials as illustrated in Figure 14 could be employed to generate the distribution 
shown in Figure 15. 

25 Figure 17 shows a plot of the resulting axial component of the gradient magnetic 
field at a series of radii in the axial direction produced by the distribution of 
magnetic material shown in Figure 15. 
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Figure 1 8 shows how the sinusoidal variation of thickness in a single unidirectional 
magnetic material can provide a sinusoidal variation of magnetization for the 
purpose of shimming. 
Theory. 

5 The design of magnets for MRI and MRS is commonly based upon the forms 

of simple geometric surfaces. The two most common of these are: a structure 
composed of parallel, finite, flat surfaces between which a uniform magnetic field is 
created as illustrated in Figure 1; and a structure with an elliptical, cylindrical 
opening within which a uniform field is produced by the surrounding distribution of 
10 material as illustrated in Figure 2. 

With reference to Figure 3, the magnetic dipolar vector field (Bd) produced by 
magnetic material can be resolved into its cylindrical polar components Bd 2> Bde and 
Bd r . In these coordinates (r, 9, z), the vector field of such an infinitesimal dipole 
has the general form: 

15 

Bd 2 (r,z) = fim (2z 2 - r 2 ) / 4n (r 2 + z 2 )* 72 (la) 
Bd e (r,z) = 0 (lb) 

and 

20 Bdr(r,z) = n m 3z r / 4n (r 2 + z 2 ) 5/2 (1 c) 

where ji is the permeability of the medium in the ROI and m is the magnetic dipole 
moment. The field due to the circulation of an electron in an atom approximates 
this very closely for distances (r 2 + z 2 ) ]n that are very much greater that the radius 
25 of the orbit (typically 10* !0 m). For the purposes of passive magnet design within 
the present context this approximation holds extremely well. 

The magnetic field from matter in bulk can be derived from equation 1. For 
the axial component of the magnetic field this is given by: 
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B 2 (r.z) = n J {M z (x,y,z) (2(z-zT - [(x-xf + (y-y 1 ) 2 ]) / <Mr - rf } dt' 



(2) 



where the underlining of individual items indicates that they are vector quantities 
and M z is the axial component of magnetization (M). Similar expressions can be 
derived for the field in the other directions. 

For paramagnetic and diamagnetic materials the magnetization M is given by: 



where Bm is the main, magnetizing field and X m is the susceptibility. In general, 
X m is a tensor of rank two. However, for isotropic materials it reduces to one of 
rank zero, e.g., for diamagnetic copper Xm = -1.0 x 10" 5 ; and for paramagnetic 
aluminium X m = 2.3 x 10" 5 . 

Equation (2) shows that the magnetic field is the convolution of the 
magnetization distribution with a spatial dispersion function. Problems of this 
nature can be conveniently solved by Fourier transforming into the kernel space (k- 
space) where the convolution becomes a simple multiplication, which greatly 
facilitates subsequent calculations. Indeed, for certain simple geometric structures, 
such as cylinders, there is no simple analytic expression for the magnetic field in 
real space, while there is one for the field in k- space. The existence of such an 
expression in k-space expedites the process of handling the above equations but it is 
not a necessary condition of the design method of the present invention. Indeed, the 
method is generally applicable, as is illustrated in the following examples. 



The design of a permanent magnetic structure according to the invention will 
now be described in general terms. The geometry and overall dimensions of the 
desired structure are set, as is a hypothetical distribution of magnetic material over 
the surface of the structure for each member of a finite set of spatial frequency 



M = Xm Bm 



(3) 



DESCRIPTION OF THE PREFERRED EMBODIMENTS 
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harmonics of unit magnetization amplitude. The field within the ROI from an 
element of the magnetic material of unit magnetization can then be determined. 
When numerical methods are employed the field must be determined to a point 
where it approximates zero in the directions in which it is to be Fourier Transformed 
5 (FT'ed). This field is then FT'ed. The resulting field in Fourier kernel space (k- 
space) is multiplied by the FT's of a set of the permissible harmonics (the lowest 
spatial frequency harmonics of the initial magnetic field distribution) of the 
magnetization on the surface. These resulting fields in k-space are subsequently 
inverse FT'ed to yield the fields in real space that are produced by each of the 

10 permissible harmonics in the selected set. These fields are fitted to the required 
field within the ROI to determine the amplitudes of the magnetization harmonics 
that minimize the sum of the squares of the deviation of the actual field from the 
required one. Fitting is performed using sets containing successively higher order 
harmonic components until the sum of the squares of the deviations of the 

15 achievable field from the desired one becomes sufficiently small. Eventually a 
point will be reached where deviations due to manufacturing tolerances will 
dominate, and no further practical improvement will be possible. 

For simple geometric structures, such as a circular cylinder, this process can 
be facilitated by directly calculating the fields in Fourier k-space. 

20 The advantage of least squares fitting lies in its convenient mathematical 

properties. For example, it commutes with the Fourier transform; as the FT itself is 
a specific example of least squares minimization. 

The following discussion serves to illustrate the method of the present 
invention only and is not intended to limit its general applicability. 

25 In a first embodiment of the present invention, the method is applied to design 

a magnet where magnetic material is distributed on the curved surface of a cylinder 
to produce a uniform field. 

For the present purposes it will be supposed that- it is required to create a 
uniform axial magnetic field B 2 of 0.1 T over a cylindrical ROI of length 0.4m and 



WO 01/29570 PCT/GB00/03905 



radius 0.125m by distributing magnetic material on the curved surface of a circular 
cylinder, such as that shown in Figure 4, having a length L of 2m and radius a of 
0.25m. The cylindrical ROI and surrounding cylindrical structure are concentric 
and coaxial along the z axis. For simplicity it will be assumed that the thickness of 
5 the magnetic shell is small compared to the diameter of the magnet, for the present 
purposes infinitesimal. In practice, a more elaborate calculation may be required to 
take into account the thickness, as described later. 

Conveniently, optimization is performed over a representative set of radii 
spanning the ROI, e.g. six radii at intervals of 20% from 0 to 100% of the radius of 
10 the ROI within the ROI. However, for the present example, optimization over just 
two curved cylindrical surfaces of diameter of 0m and 0.25 m and lengths 0.4m, was 
sufficient, as demonstrated by table 3 and Figure 12 (explained below). 

The distribution of magnetic material on the surface of the structure shown 
in figure 4 can be conveniently described in terms of a series of orthogonal 
15 azimuthal and axial, sinusoidal and cosinusoidal, harmonics. The azimuthal 
harmonics form a discrete set of cosinusoidal terms cos (m<j> + a) and the axial 
harmonics form a truncated set of H(z) cos (n9 + p), as illustrated in Figure 5, 
where H(z) is the unit step function given by: 

20 for |z| < L/2 H(z) = 1 

(4) 

for |z| > U2 H(z) = 0 

where L is the length of the cylindrical surface and a and P are phase terms that take 
25 account of both the sinusoidal and cosinusoidal terms. 

In many applications the set of harmonics is severely curtailed by 

considerations of symmetry. For example, for the present uniform field there is no 

azimuthal variation and all azimuthal terms other than m=0 are zero. 

Consequently, the problem of optimization reduces to a one-dimensional one of 
30 determining just the axial variation. Further, as the ROI and the structure are 
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concentric the resulting even symmetry of the system implies that all sinusoidal 
terms are also zero leaving only the purely cosinusoidal terms. It remains to 
determine the amplitudes of each of the significant cosinusoidal terms to produce an 
optimum uniformity of the field within the ROI on the basis of the minimum sum of 
5 least squares deviations of the field. 

Now the field from an elemental magnetic dipolar ring of unit magnetization 
can be determined by numerically integrating equation (2). Figure 6 shows how the 
resulting magnetic field varies with distance, in the axial direction, at a radius of 0.1 
m from such a dipolar ring of diameter 0.5 m. 
10 From equation (2) we also know that the magnetic field from individual axial 

harmonics is the convolution of the field from this ring with the harmonic. This 
convolution is most conveniently performed by simple multiplication in the Fourier 
domain. Consequently, the next step in the design process is to FT this field in the 
axial direction z; the resulting FT is shown in Figure 7. 
15 Alternatively, for a simple geometric surface such as this, the FT in the axial 

direction of the axial component Br z of the ring field, Br, can be expressed 
analytically as: 

forr<a Br z (r,m,k) = - u a k 2 I m (|k| r) K m (|k|a) / 2n (5a) 

20 forr>a Br z (r,m,k) = - ua k 2 I m (|k| a) K m (|k|r) / 2n (5b) 

where I m and K m are the modified Bessel functions. 

Figure 8 shows the FT in the axial direction of the axial component Br z of the 
ring field, Br, evaluated directly by equation 5a. This is substantially identical to the 
25 field derived by numerically integrating the generally applicable equation (2) and 
then Fourier Transforming it, as can be seen by comparing Figures 7 and 8. 

Although evaluating the FT directly by means of equation (5a) facilitates the 
process, this method is not applicable to structures where an exact solution in the 
Fourier kernel domain is not available. The longer route of numerical integrating 
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equation (2) and subsequent FT'ing has general application for any structure, as it 
does not rely upon there being an analytic form of the equations in either real or k- 
space. 

Now, for an open cylinder of restricted length, L, the total magnetization, M z , 
5 can be expressed as the sum of the permissible axial harmonic modes, n, on the 
surface as follows: 

M z (a, z) = I„ C n H(z) cos (2 it n z / L - + ) (6) 

10 where C„ is the amplitude of the nth magnetization harmonic. The FT of this in the 
z-direction is: 

M z (a, k) = 7i L Z n C n [e* sinc(k U2 - mr) + e'* sinc(k L/2 + mi)] (7) 



15 



25 



is 



Therefore, the FT of the total resultant field from the cylindrical surface i 
simply the product of the FT of the magnetization with the FT of the elementary 
ring field within this domain, i.e.: 



20 B z(r. m, k) = M 2 (a, k) Br 2 (r, m, k) 



(8) 



Optimization of the distribution can is then solved by least squares 
minimization of the first few terms (C n ) over the ROI in real, as against kernel, 
space. The number of terms is dependent upon the geometries of the coil structure 
and the ROI. This can be determined by solving with increasing numbers of terms 
(C„) until a sufficiently homogeneous field has been attained. 

In the above example the first 3, 4 and 5 harmonic amplitudes (coefficients) 
were employed and the worst deviation from the field in the ROI determined (see 
Table 1). Optimization over more than 5 harmonic coefficients has no further 
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significant effect As can be seen from Table 1, the maximum deviation of the field 
in the ROI drops to ~ lppm within the first 5 terms. In any case, practical 
considerations, such as manufacturing tolerances and the finite thickness of the 
structure, mean that further improvements in field uniformity are impractical. 
5 Further approximation can be achieved by shimming the magnet in some fashion, 
for example along the lines discussed in the second embodiment below. 



Table 1 : Field deviation with number of harmonic components employed. 



No. of Harmonics 


Range of Harmonics 


Max deviation in ROI 


3 


0-2 


200 ppm 


4 


0-3 


16 ppm 


5 


0-4 


1 ppm 



10 The harmonic amplitudes of table 2 were used to determine the axial 

distribution of magnetic material (M = 1.28 x 10 6 amp/m, equivalent to the 
magnetization of saturated iron) on the surface of the hollow cylindrical structure, as 
shown in Figure 9. 

Figure 10 shows a practical implementation of a magnetic structure designed 
15 in accordance with the invention, to generate the distribution shown in Figure 9. 



Table 2: The harmonic amplitudes for a uniform field. 



Harmonic No. 


Amplitude (amps/m) 


0 


-1,143,966 


1 


-220,943 


2 


37,362 


3 


-6,821 


4 


749 
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Figure 11 shows the resulting magnetic field at a series of radii (0, 2.5, 5, 7.5, 
10 and 12.5 cm) plotted in the axial direction produced by the distribution of 
magnetic material, according to the present invention, as shown in Figure 9. Figure 
12 shows the deviation of the above magnetic field from the desired, uniform one, 
5 in parts per million (ppm), over the ROL Figure 13 shows how the field diverges 
rapidly outside the ROI in the axial direction, but remains reasonably homogeneous 
in the radial direction; see also Table 3 below. 



Table 3: Maximum field deviation in the ROI at each radius. 



Radius cm (% fraction of magnet) 


Maximum deviation in ROI 


0.0 ( 0) 


0.460 ppm 


2.5 (10) 


0.160 ppm 


5.0 (20) 


-0.541 ppm 


7.5 (30) 


-L133 ppm 


10.0 (40) 


-1.011 ppm 


12.5 (50) 


0. 1 64 ppm 


15.0 (60) 


2.011 ppm 


20.0 (70) 


-4.730 ppm 



10 

As demonstrated by the above figures, the resulting field is uniform to a very high 
degree in the central region. 

Figure 1 1 serves to highlight the fundamental inapplicability of methods that 
attempt to invert the field equations to the design of permanent magnetic structures. 

15 Due to the complexity of the field, it is inconceivable that even someone highly 
conversant with the design of magnets could arbitrarily specify a field such as that 
of Figure 1 1 with the required degree of accuracy to ensure a practically realizable 
distribution that would yield a reasonably good approximation to the desired field in 
the ROI; essentially this would require a knowledge of the result of the process 

20 before the process had begun. Certainly, this would not yield the uniformity attained 
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by the present technique which only requires the specification of the known, desired 
field within the ROl. 

In accordance with a second embodiment of the present invention, the design 
method is applied to homogenize an existing field by distributing magnetic material 
5 on the curved surface of a circular cylinder, of the type shown in Figure 4. For the 
purpose of illustrating this example, the existing field is assumed to be oriented in 
the z-direction and is substantially uniform except for a first order (linear) 
transverse variation (inhomogeneity) in the ROL It is to be appreciated that the 
design method is by no means limited to compensating for a linear variation, but is 
10 generally applicable to shimming for any inhomogeneity. 

As in the previous embodiment, the magnetic surface is arbitrarily assumed to 
be 0.5 m in diameter and 2 m in length, with a cylindrical ROI of diameter 0.25 m, 
and length 0.4 m. Optimization was performed on a representative curved 
cylindrical surface of diameter 0.2 m and length 0.4m to produce a transverse linear 
15 gradient field of 0.1 milliTesla per metre. Again, table 5 below and Figure 17 
demonstrate that this is indeed representative of the ROI. 

For simplicity it will again be assumed that the thickness of the magnetic shell 
is small compared to the diameter of the magnet, for the present purposes 
infinitesimal. In practice, for the creation of shim structures this approximation will 
20 hold and more elaborate calculations will generally not be necessary. 

From symmetry considerations, an azimuthal variation of cos (<f>), 
corresponding to a linear transverse variation, in the field can only be produced by 
a first order variation (m=l) in the magnetic material. This again reduces the 
problem to a one-dimensional one of determining just the axial variation. 
25 The FT of the field can be determined directly from equation 5a. Optimization 

of the shim is then conveniently solved by least squares minimization of the first 
few terms (GO over the ROI in real, as against kernel, space. 

In this embodiment coefficients 1, 2, 3 and 4 were necessary as can be seen 
from table 4. From inspection, it was noted that the zeroth coefficient contributed 
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negligibly to the central field, so this too could be omitted. The amplitude of the 
required axial distribution is shown in Figure 15. The distribution varies azimuthally 
as cos(0) (with amplitude = 1.28 x 10 6 amp/m) on the surface of the cylindrical 
structure to provide a magnetic field that varies optimally, linearly, transversely 
5 (Bz(x) = x. 1 00 microTesla/m) over the ROI. 

Figure 16 shows how a shim magnetic structure designed in accordance with 
the invention and using diamagnetic and paramagnetic materials as illustrated in 
Figure 14 may be implemented in practice, to generate the distribution shown in 
Figure 15. 

10 Figure 17 shows the resulting gradient magnetic field at a series of radii (0, 

2.5, 5, 7.5, 10 and 12.5 cm, 15 cm) , plotted in the axial direction, produced by the 
distribution of magnetic material shown in Figure 15. Note that the linearity remains 
excellent far outside the ROI. 

15 Table 4: The harmonic amplitudes for a linear transverse field. 



Harmonic No. 


Amplitude (amps/m) 


1 


247 


2 


-59 


3 


15 


4 


-2 



The resulting magnetic field is linear to a very high degree in the central 
region, see table 5 and figure 17. Certainly, the degree of linearity is well within the 
limitations of first order approximation adopted and of the errors associated with 
20 manufacturing. Consequently, such a shim designed in accordance with the present 
invention will correct the linear inhomogeneity to a very high order of accuracy. As 
mentioned above, the design method of the present invention is not limited to 
correction of linear variations, and a similar to degree of correction can be 
anticipated generally, for any conceivable inhomogeneity. 
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Table 5: Center field in the ROI at each radius for a linear transverse shim. 



Radius cm (% fraction of magnet) 


Bz (microTesla) 


0.0 (0) 


0.000000 


2.5 (10) 


2.499998 


5.0 (20) 


4.999995 


7.5 (30) 


7.499989 


10.0 (40) 


9.999974 


12.5 (50) 


12.499929 



It is to be appreciated that, generally, the lowest spatial frequency harmonics 
5 are chosen, but this will not always be necessary as witnessed by the example 
above. 

The number of terms over which the optimization is performed may be 
reduced by incorporation of an additional, hypothetical, arbitrary distribution of 
magnetic material over the surface of the structure into the set of spatial frequency 
10 harmonics. By "arbitrary" it is to be understood that a suitable distribution is 
determined by appropriate consideration of the desired field. For example, to 
achieve a field that varies linearly along a flat surface of magnetic material, an 
appropriate distribution might be one whose magnetization also varies linearly 
along that surface. 

15 The method of the invention can be applied to the design of shims in several 

ways. One of these is to fit the dipole harmonics directly to the entire field 
inhomogeneities, as illustrated above. Another is to fit the dipole harmonics to 
individual spherical harmonic terms of the field. The rationale for the latter is that it 
is common practice to map the field in these linearly independent spherical 

20 harmonic terms. The required shim is then the linear combinations of the shims for 
the individual spherical harmonic terms weighted appropriately. This second 
method is advantageous when shimming a series of magnets of the same physical 



- 15- 



WO 01/29570 



PCT/GBOO/03905 



dimensions, as the shim dimensions remain constant and consequently the spherical 
harmonic shims form a standard set. The resulting shim from the example above 
effectively also represents one component of such a set. 

The required distribution of magnetization can be achieved in several ways. In 
5 one preferred embodiment the thickness of a single magnetic material can be varied 
in proportion to the required distribution. This is perfectly adequate for the 
production of uniform fields such as described in the first embodiment provided the 
thickness of the material is small in comparison with the sizes of the structure and 
the ROL Where this is not the case the process must be refined by shimming along 

10 the lines outlined for the second embodiment. 

For the purpose of shimming the use of a single unidirectional material cannot 
produce a magnetization that varies sinusoidally, but if an additional uniform field 
can be tolerated simply varying the thickness of the material sinusoidally in the 
manner shown in Figure 18 produces acceptable permissible harmonic 

15 magnetization components of the form A-cos(mG) and A-sin(mG), in the example 
illustrated, 2-sin(30). 

In another preferred embodiment distributing materials of differing 
magnetization in the vicinity of the ROI can be used to achieve the desired local 
magnetization. By using differing thicknesses of materials that have opposite 

20 susceptibilities harmonic components of the form cos(mG) and sin(m0) are possible. 

For example, for a cylindrical surface consisting of a semicircular (0 < 0 < n) 
distribution of paramagnetic aluminium (xm = 2.3x1 0" 5 ), whose thickness varies as 
0.435 A |sin(0)| and an opposing semicircular (-7i<0<O) distribution of 
diamagnetic copper (Xm = -1.0xl0~ 5 ) whose thickness varies as A |sin(0)|, produces a 

25 magnetization that varies as A sin(0) to first order over the whole surface 
(-71 < 0 <ti), in a similar manner to that shown in Figure 14. Alternatively, the 
magnetization distribution shown in Figure 1 5 could be achieved by machining the 
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above sinusoidal contour on cylinders of paramagnetic and diamagnetic material in 
this manner. 

A material of any required magnetism can be constructed by combining 
paramagnetic and diamagnetic materials in the required proportions. For example, a 
5 material that has effectively zero magnetization can be formed from a combination 
of 30.3% aluminium and 69.7% copper. This could be a great value in 
manufacturing objects for use in magnetic fields that would not affect the field, or 
be affected by the field. In practice, the proportions may have to be altered slightly 
to take into account the higher order effects when such materials are combined. 
10 For either creating a uniform field or altering an existing one, the 

thickness of the magnet will be a significant fraction of the diameter and the 
consequent placement of material will affect the actual attainable field. To take 
account of this, the method can be applied iteratively as follows, by: 

-calculating the distribution of magnetic material to produce an optimum 
15 approximation to a desired magnetic field in the ROI, in accordance with the 
method described herein 

-determining the total resultant field in the ROI taking into account the 
thicknesses and placements of magnetic materials involved 

-calculating a corrected magnetic material distribution representative of the 
20 difference between said total resultant field and said desired field 

-repeating the latter two steps until the field attained approximates the 
desired one in the ROI to a sufficient degree. 

It is to be appreciated that the homogeneity of the field generated by a 
magnetic structure designed in accordance with the present invention can be further 
25 enhanced by additional active shimming. Resistive electromagnetic shimming may 
subsequently be employed to further improve the homogeneity during the use of the 
magnet. This is particularly relevant at high field strengths with gradient echo 
techniques where the natural susceptibility of the sample distorts the field, in other 
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words, where resistive electromagnetic shimming can be employed to improve the 
homogeneity for each sample. 

Whilst the above embodiments have been described in the context of magnets 
of cylindrical form, it is appreciated that the method of the invention is equally 
5 applicable to the design of magnets of other geometric forms, one example of which 
is the parallel flat plate arrangement shown in Figure 1 . 
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